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ABSTRACT 

Consistency of present-day lattice QCD simulations with dynamical 
( "sea" ) staggered fermions requires that the determinant of the staggered- 
fermion Dirac operator, det(i5), be equal to det^(-Drg) det(T) where 
is a local one-flavor lattice Dirac operator, and T is a local operator con- 
taining only excitations with masses of the order of the cutoff. Using 
renormalization-group (RG) block transformations I show that, in the 
limit of infinitely many RG steps, the required decomposition exists for 
the free staggered operator in the "flavor representation." The resulting 
one-flavor Dirac operator D^-g satisfies the Ginsparg- Wilson relation in the 
massless case. I discuss the generalization of this result to the interacting 
theory. 



1. Introduction 

Lattice QCD simulations with dynamical ( "sea" ) staggered fermions PP are pro- 
viding predictions for liadronic observables with unprecedented accuracy ^ . In these 
numerical calculations, it is crucial that all sources of error be under systematic con- 
trol. This raises the question of the validity of the "fourth root trick" used in these 
simulations. 

Let me briefly explain the problem (more details may be found e.g. in ref. In 
four dimensions, the staggered Dirac operator D is a one-component lattice operator 
which, in the free-field case, has sixteen poles in the Brillouin zone. These poles 
combine into four Dirac fermions (with a total of sixteen degrees of freedom) in the 
continuum limit. To account for three dynamical quarks - up, down and strange - 
the Boltzmann weight used for generating the dynamical configurations involves the 
factor^ 

det^/'(D(m„)) det^^\D{md)) det^/'(D(m,)) . (1) 

Taking the fourth root of each staggered-fermion determinant ensures that the lattice 
theory describes three (and not twelve) quarks in the continuum limit. While the 
"fourth root trick" is necessary in practice in order to reach the desired continuum 
theory, it is not obvious that this trick is consistent. The question is whether the 
gauge-field configurations generated with this Boltzmann weight correspond to a local 
lattice theory. If one could show that det{D) is equal to the fourth power of the 
determinant of some local one-flavor lattice Dirac operator,^ this would provide a 
positive answer. 

In fact, a far weaker condition is sufficient to guarantee locality, and, hence, 
consistency.^ Suppose one can show that 

det{D) = det^(Ag) det(r) . (2) 

Here Dj-g is a local one-flavor lattice Dirac operator, and T is a local operator con- 
taining only excitations with masses of the order of the cutoff. We may now write 
det(T) = exp(— S'efr(f/)). Since T contains only excitations with cutoff masses, we ex- 
pect that the effective action Ses{U) is local; trivially, the same is true for (l/4)S'eflf(t/). 
By eq. the fourth-root trick then amounts to using dynamical fermions with the 
local Dirac operator D^-g, together with the modification of the gauge-field action by 

(l/4)5eff(f/).^ 

The natural framework to realize relation Q is through RG block transforma- 
tions. After introducing the relevant concepts in Sect. 2, the free staggered-fermion 

"'^In practice, m„ = in the simulations. Note also that the simulations correspond to a 
"hybrid" and/or partially-quenched theory, since the "sea" and the "valence" quarks may differ in 
their masses, as well as in the details of the discretization used for each. 

^In this paper, a local operator means an operator whose kernel decays exponentially with the 
separation |a; — ?/|, with a decay rate which is 0(1) in lattice units. A similar notion of locality 
applies to the effective action Ses{U) discussed below Eq. 

■^This observation was recently made in ref. 

^Locality of det^^^(T) in the free theory is addressed in Sect. 3. 
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operator is dealt with in Sect. 3, which constitutes the main part of this work. Using 
the "flavor representation" of staggered fermions |5j it is shown that the decompo- 
sition (j21) is realized in the limit of infinitely many RG blocking steps. Central to 
this discussion is a theorem on the locality of RG-blocked Wilson fermions proved 
in ref. 0; only a trivial amendment is needed in order to generalize the theorem to 
the flavor representation of staggered fermions. The limiting operator is constructed 
explicitly. (For Wilson fermions, see ref. [7^. See also ref. |Bj for an RG treatment of 
staggered fermions within the one-component formalism.) Like any fixed-point oper- 
ator [9;, the limiting operator satisfies the Ginsparg- Wilson (GW) relation [TUl lllj . 
Sect. 4 contains a discussion of some of the issues (both theoretical and practical) 
that arise in the interacting theory. 



2. Renormalizat ion-group transformations 

We first consider a general setup, following refs. [^IHl- Starting from a bilinear 
fermion action with Dirac operator Dq, an RG blocking is introduced via the identities 



d^^jJd^p exp{—ipDoijj) (3a) 
dipdipdxdx exp {-ipD^ip - a{x - ^<5"^)(x - W)) (3b) 
dipdijjdxdxd'qdf] exp [—^pDoip + VV + ix ~ ^Q^)v + vix ~ Q'^))(3c) 
= det-\Gi) J dxdx expi-xD^x) ■ (3d) 

Here live on the "fine" lattice with spacing oq, whereas XtX (and the auxiliary 
field Tj, rf) live on the "coarse" lattice with spacing ai equal to an integer multiple of 
oq. In this paper we set ai = 2aQ. The RG blocking kernel Q is a rectangular matrix 
satisfying 

QQ^ = cl , (4) 
where / is the identity matrix on the coarse lattice. Explicitly, 

Gi = T, = {Do + aQ^Q)-\ (5a) 
Di = a-a^QGiQ\ (5b) 
Di' = a-^ + QD^'QK (5c) 



Eqs. fl5a|) and (|5b|) are obtained by integrating over ip,^ in Eq. (j3b|) . and Eq. (jog) 

by first integrating over ipjip and then over 77,77 in Eq. (jHcj) . Iterating the blocking 
transformation we have 

r,- = (D,_i + ag^-'^tgO))-! , (6a) 

Dj = a - a2g(^)rjg(j)t , (6b) 

Dj' = + g(^)D-_\g(^)t^ (6c) 
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where Q^^^ denotes the blocking kernel from the (j — l)-th lattice (with spacing 2^~^aQ) 
to the j-th lattice (with spacing 2-'ao). We may also go directly from the finest to the 
coarsest lattice: 

Gn = {Do + anQiQn)~\ (7a) 

1 — c 

Dn = an- alQnGnQl , = «- , (7b) 

1 — c" 

= a-' + QnD^'Qi, (7c) 

where = Q^'^^Q^'^~^^ ■ ■ ■ Q^^\ Eq. (f7c|) follows by iterating Eq. (jUc|) while using 
Eq. (0)). The other equations follow by noting that the product of n blocking trans- 
formations can also be represented as a single "big" blocking transformation as in 
Eq. (0), provided we let x^X live on the coarsest lattice, and we make the replace- 
ments a an, Q Qn- Hence each relation in Eq. must match the corresponding 
one in Eq. 



3. Free staggered fermions 

We now turn to free staggered fermions. In the flavor representation [5^ the 
staggered Dirac operator has four-component spin and flavor indices, and is given 
explicitly by 

Do = a-^ J2 ((^M ® + (75 ® T-5r^)A^) + m , (8) 

where a is the lattice spacing, V^/(a:) = {f{x + fi) — f{x — fi))/2, and A^/(x) = 
(2/(a;) — f{x + /t) — f{x — fi))/2. The usual Dirac matrices are 7^, while the 
constitute another set of Dirac matrices acting on the flavor index. (Taken together, 
(7^ J) and (75 ® T^) form a representation of the eight-dimensional Dirac algebra.) 
For m = 0, Dq is anti-hermitian. Going to momentum space one has L>o ^ = ^'^dI 
where 

Do = a~-^ ((7^ I)i sin(p^a) + (75 nT^){l - cos(p^a)) + m , 
n = dIDo = a-^^{ sin^(p^a) + (1 - cos(p^a))^) + . (9) 

We will apply n block transformations to the Dirac operator (jSI). We set m = 0, 
hence we may write 

L>o"' = D^\p- a) = -J2 (^(7m ® I)Al{p; a) + (75 ® nT,)Bl{p; a)) . (10) 

We will hold fixed the lattice spacing obtained in the n-th step. We thus set 2"ao = 1, 
or ao = 2~". The blocking kernel Q'--'^ is defined as follows. We label the sites of the 
(j — l)-th lattice by four integers I = (/i, Z2, ^3, ^4), & Z. A site /' = (Z'^^, /g, ^3, ^4) 
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of the j-th lattice is identified with the site 21' = (2/1,2/2,2/3,2/4) on the {j — l)-th 
lattice. The blocking transformation assigns to a field variable on the j-th lattice its 
arithmetic mean over a 2^ hypercube on the (j — l)-th lattice. Explicitly, {Qf){l') = 
Q(/', /)/(/) = 2-^T.r^=o,i /(2^' + r). This definition implies c = 2"^ in Eq. Q. 
Using eq. (fTcj) we obtain 

D-\p) - 
Kip) - 

b;{p) - 

IQn(p,A:("))P = 

where -vr < < vr and for each /i, k],""' = -T~^,-T-^ + l,...,2""i - 1. To 
arrive at Eq. ()lld|) . observe that we have set a = 1 for the coarse-lattice spacing, 
so that a single-step blocking kernel is a mapping into this lattice from a lattice 
with spacing equal to 1/2. For each /z, this kernel has the momentum representation 
(l/2)(exp(zp^/2) + 1) = exp(ip^/4) cos(p^/4) = exp(ip^/4) sin(p^/2)/(2 sin(p^/4)). 
Eq. (jlldj) follows by going "backwards" down to the fine lattice with spacing ao = 2^". 

The massless staggered operator satisfies {-Do, (75 ® ^5)} = 0, and therefore the 
staggered action is invariant under a t/(l) chiral symmetry.^ For the RG-blocked 
operator we have Eq. pia|) . which implies that D~^(a;, y) anti-commutes with (75(^X5) 
except for x = y. In fact, 

{D„(75®r5)} = 2a^iZ}„(75®r5)D„, n>l. (12) 

This is recognized as a generalization of the GW relation ^U]. Thus, after the 
first blocking transformation the U(l) chiral symmetry gets modified to a Ginsparg- 
Wilson-Liischer (GWL) chiral symmetry [TT] . 

For Wilson fermions, it was proved in ref. |E] that the RG-blocked operator 

Dr,{x-y)= [ dpe'P^^-y^D^{p), (13) 

J BZ 

is local. Here Jg^ dp = J^^ denotes the integration over the Brillouin zone of the 
coarse lattice. As explained earlier, this means that Dn{x — y) decays exponentially 
with \x — y\, and the decay rate is 0(1) in units of the coarse-lattice spacing. The 
bounds established in the course of the proof are uniform in n, and hold for a < a 

^This symmetry correspond to the J7(l)e symmetry of the one-component formahsm [ISj. Note 
that the interpretation of this symmetry in the continuum hmit - axial, vector or some combination 
of them - depends in general on the choice of the staggered mass term For the simple mass 
term of Eq. it is a chiral symmetry. 



- E (^(^^ ® ^^-^W + ® r,T,)B;{p)) , (11a) 

J2 + 27rA;("); 2-")|g„(p, A;("))p , (lib) 

J2 KiP + 27rA;("); 2-")|Q„(p, fc("))p , (11c) 

sin(p,,/2) 



n 



2" sin((p^ + 27rA;^,"V2''+^ 



lid) 
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where d > is an 0(1) constant whose actual value can be worked out by keeping 
track of the details of the proof. 

We now argue that the proof of locality continues to hold if, at the starting point, 
we replace the Wilson operator by the staggered Dirac operator (jH)). This amounts 
to replacing the Wilson term W hj a "skewed" Wilson term Wgt, where 

W = ^(1 - cos(p^a)) , 
W,t = 5]](75®r5r^)(l-cos(p^a)). (14) 

The proof requires lower and upper bounds on as a function of p^a. Introducing the 
vector-space norms I I ^ = ( jx^l''') we observe that = \W\ = |1— cos(p^a)|i. 
In the staggered case, we have the operator norm || Wst |P = VTstW^jt = ~ 
cos(p^a))^, or equivalently || Wst || = |1 — cos(p^a)|2. Since the following equivalence- 
of-norms inequalities hold in d dimensions 

rf"'/'|xj2 < \x^\i<d'/^\x^\2, (15) 

it follows that every lower or upper bound on \W\ entails a corresponding bound 
on II Wst II; and vice versa. This simple argument shows that, indeed, the proof 
given in ref. [H] generalizes to the RG-blocked staggered Dirac operator in its flavor 
representation. 

The key physical input that goes into the proof [Hj is that Dq^{p) and D^^{p) 
share the same singularity as p — 0, namely, —i^^l^Vn/v^- Indeed, the singularity 
of D~^{p) arises only from the /c^"'' = term on the right-hand side of Eq. ()llb|) . 
Factoring out this singularity by writing -Do(p) = Rn{p)Dn{p), one can prove that the 
operator Rnip) is analytic in p^ and that both Rn{p) and R~^{p) are bounded. This 
is then used to prove that Dn{p) and Gn{p) are analytic, and that D^ip), Gn{p) and 
G~^{p) are bounded. Exponential localization of the corresponding coordinate-space 
kernels follows from general theorems. 

We will next show that, in the limit of infinitely many RG steps, D„ becomes 
diagonal in flavor space. The flavor-mixing part of D~^{p) is given by Eq. (jllc|) . 
where explicitly 

B'ip + 2vrA:W- 2-) = 2 sin^((p, + 2vrA:("))/2»+^) 

^ ' 2«^^ {sm\{p, + 2nkjr^)/2-)+Asm\{p, + 2nki''^)/2-+^)) 

(16) 

It is easy to see that |i3°(p + 27ik^"-'>; 2 ")| < Ci2 " where Ci = 0(1). Away from the 
singularity at p + 27ik^"''^ = this is evident. For \p + 27rfc*^")| -C 2", the same result 
follows using sin(x) ~ x for X < 1. Thus, |i3°(p -t- 27rfc("); 2-")| = 0(2""), uniformly 

in p and k^'^K In addition, \Qn{p,k^"'^)\'^ < C2Yl^\pu\'^ /{{pu + 27rA;^"-*p where again 
C2 = 0(1) jH]. Hence the A;^"''-summation converges, and |i3JJ(p)| = 0(2"") for all p. 
Taking the limit n — > cxd we conclude that B'^{p) 0, uniformly in p. 



6 



The inverse fixed-point operator obtained in the hmit n — >■ oo can be expressed 

D^\p) = D;^\p) (^I=(^a^'-^J2 l.A^ip^j ® / , (17) 

where Ooo = (15/16)a, and A'^{p) ~ P^/p^ for <^ 1. This shows that D^o is 
diagonal in flavor space and satisfies the (standard) GW relation.^ Letting Goo = 
lim^^ooGn, Eqs. (jSI) and ^7i) imply 

lim det^/^(Do(ao = 2"")) = det(Ag) det^/^(G^^) . (18) 

n^oo 

In view of the locality and boundedness properties estabhshed above, the desired 
decomposition is achieved in this limit! 

The essence of the RG blocking is that it distills the long-distance dynamics, 
extracting it out of the underlying short-distance theory. The long-distance dynamics 
is contained in Dj-g which is manifestly diagonal in flavor space, while all the flavor- 
mixing effects are contained in G^. Since is analytic in momentum space and 
has an 0(1) gap, its fourth root shares similar properties. Hence Gn^^^ is local, and 
has only cutoff-mass excitations, uniformly in n. 

Let us elaborate on this last statement. Observing that a has mass dimension 
equal to one, and focusing e.g. on the first blocking step, an RG transformation 
works by first replacing Dq with Gi^ = Dq + aQ''Q {cf. Eq. (jSEJ). Now, the massless 
operator Dq has vanishingly small eigenvalues near = 0. The contribution from the 
blocking kernel, aQ^Q, lifts these small eigenvalues and generates an 0{a) = 0(1) 
gap in the spectrum of G^^. We may define the fourth- root operator for any finite n 
via 

Mnix,y) = / dpe'P^^-y^MniP) , MniP) = (G„(p)Ot (p))-i/8 _ ^^g^ 
Jbz 

In this equation, x, y take values on the fine lattice, and jj^^ dp denotes the integration 
over the fine-lattice Brillouin zone. The argument why Ain{x,y) is local is standard 
jH]. If we let one of the momentum components become complex, the singularity 
closest to the real axis will be at a distance which is 0(a„). Deforming the contour 
of integration, this implies that A4n{x,y) decays exponentially with \x — y\, with a 
decay rate which is 0(a„), namely, 0(1) in units of the coarse-lattice spacing. 

It is interesting to compare this result to ref. j3] which attempts to find an operator 
jV such that det(A/') = det^'''*(Do); without the help of RG transformations (see also 
ref. [in]). In this case, the gap is provided by the physical mass, and the decay rate of 
the square-root kernel needed in the construction is found to be mj a). Hence, 
the limit m — > is problematic. 

With the help of RG blocking, the small-distance scale relevant for the fourth 
root is ao, the spacing of the original fine lattice, and the relevant large- distance scale 

^The relevance of the GW relation for establishing consistency of the fourth- root trick was recently 
pointed out in ref. |15| . 
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is a, the spacing of the coarse lattice. In comparison, in refs. fHl H] the relevant 
short- and long-distances scales are a and 1/m respectively. The use of RG blocking 
effectively achieves the replacements a —>■ a^, l/m —>■ a. While the mathematics is 
similar, the physical conclusion is different. The "fourth-root" kernel Ain{S:,y) is 
long-ranged with respect to the fine- lattice scale Oq (in analogy with refs. jlHlllI); but 
the same kernel is short-ranged with respect to the coarse-lattice scale. This is true 
uniformly in n, hence also in the limit n —>■ oo. 

The RG analysis generalizes to the case that the mass term is switched back on 
in Eq. ^7 We find that Eq. ()lla|) still holds, except that the explicit expressions for 
A^{p), B^{p) and get modified (now = a^^(p) becomes a non-trivial function 
of p). The proof of locality generalizes to m 7^ for Wilson fermions 0, and the 
same should be true for the flavor-representation staggered fermions. Also, clearly 
B^{p] m 7^ 0) < i3°(p; m = 0). Therefore, for — > 00, B^{p) tends to zero as before. 
The limiting operator Dj-g, defined by the first equality in Eq. (fTTj). is again diagonal 
in flavor space. Of course, for m 7^ 0, D^g will not satisfy the GW relation any more. 

4. Interacting staggered fermions 

It is unlikely that rigorous theorems such as those of ref. (H] will ever be gener- 
alized to an interacting lattice theory. Still, physical intuition suggests that similar 
statements on locality and boundedness may hold true in an interacting theory as 
well. In this section, I address some of the issues that arise when dealing with an 
interacting theory. 

The question is what are the properties of an RG-blocked lattice theory, when the 
initial interacting theory involves one-component staggered fermions, and the fourth- 
root trick is applied. In an interacting theory RG transformations may be realized in 
numerous ways. A common feature is that RG transformations naturally give rise to 
a "two-cutoff" theory: the RG-blocked theory living on a coarse lattice with spacing 
a is obtained after applying n RG transformations to an initial theory defined on a 
fine lattice with spacing -C a. For example, in the context of "perfect action" one 
applies RG transformations to fermion and gauge fields alike, and the limit n — >■ 00 
(and Oo — ^ 0) is taken while keeping a fixed |9|. 

Here I will limit the discussion to a simpler framework, where RG transformations 
are applied only to the fermion variables. Among other things, this has the advantage 
that some simple tests can be carried out on existing dynamical configurations. 

The first problem that must be tackled is that, as discussed below, the interacting 
theory is defined using one-component staggered fermions for a reason El ■ In the 
free-field case, there is a unitary operator Qq that maps the one-component staggered 
operator, defined on a lattice with spacing oq, to the flavor-representation staggered 
operator on a lattice with spacing oi = 2ao |5J. In the interacting case, the mapping 
must preserve gauge invariance but there is no unique, obvious way to define it. 

^It is desirable (though not a necessary condition for ma ^ 1) to choose the same relative sign 
for m and a, such that m + aQ^Q is a strictly positive operator. 
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We propose to deal with this problem by a single RG blocking step which keeps 
the number of fermionic degrees of freedom unchanged. Specifically, in Eq. (jS)) we 
take to be single-component fields, while x^X ^^re four flavors (or "tastes") of 
Dirac fields. Dq is now a covariant, interacting one-component staggered operator. 
For the blocking kernel we choose some covariant version of Qq, denoted Q^^-* below, 
defined in terms of the link variables on the fine lattice. Eqs. and (jSj) provide 
explicit expressions for Di, which is the resulting staggered operator in the flavor 
representation, as well as for Gi, whose (inverse) determinant is picked up when 
performing this non-trivial change of variables. (We will shortly return to the role of 
det{Gi^).) Using the natural embedding x ^ x from the coarse lattice to the fine 
lattice (as described above Eq. (fTTj) ). both Di and Gi are gauge-covariant functions 
of the link variables on the fine lattice. 

Having thus constructed a "flavor representation" in the interacting theory, we 
may apply n additional, ordinary RG blocking transformations to the fermions. These 
block transformation dilute the number of fermionic degrees of freedom while main- 
taining gauge invariance, provided that we keep choosing blocking kernels which are 
covariant with respect to gauge transformations on the original fine lattice. This is 
naturally realized if, for any point x' and any point y such that Q^^\x',y) 7^ in 
the free theory, we construct the covariant kernel Q^^\x', y) by summing over Wilson 
lines that go from x' to y on the original fine lattice. 

An important question is what is the fate of the global symmetries of the orig- 
inal, interacting one-component staggered-fermion theory. As for the chiral U{1) 
symmetry, we have seen that it becomes a GWL symmetry (Eq. ()12|l ). Now, it 
may be impossible to preserve manifest hypercubic invariance in the construction of 
and the same is true for the staggered shift symmetry [21- Since the original 
one-component theory has exact hypercubic and shift symmetries, by Eq. ()3d|) . any 
breaking of these symmetries induced in the above-defined flavor representation by 
Dn must be exactly compensated by det(G~^). In other words, the effective action 

= log(det(Gn)) should automatically contain the local "counter-terms" needed to 
restore exact invariance. 

It is known that the attempt to construct an interacting theory directly in the 
flavor representation gives rise to induced cutoff-scale masses which, on top of that, 
violate Lorentz invariance [d- The interacting staggered theory is defined in the 
one-component formalism because its global symmetries, including in particular shift 
symmetry, forbid these disastrous mass terms ^3]. It is therefore crucial that the 
cancellation mechanism proposed above will indeed be operative when the flavor 
representation is constructed though "RG blocking" from the one-component theory. 
Interestingly, there exists a general RG-blocking result which provides direct evidence 
that the above cancellation mechanism indeed works as it should. Introduce the 
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"telescopic sum" [H] 



j=0 

n-l 

= SnD^-^Sn + ''Y^ SjTj-^^iSj , (20) 
j=0 

where we have set Qo = I and where 

S, = D,'Q]D^, S, = DjQ,D,\ (21) 

and 

f,+i = - DT^Q^^+'^^D,+,Q^^+'^DJ' = r,+i . (22) 

The equahty Fj+i = Fj+i follows by substituting Eq. (jUcj) for -Dj+i and re-expanding 
as a geometric series. 

Consider now an interacting staggered-fermion theory with light or massless 
quarks. Let us examine the long-distance behavior of the various kernels in Eq. ()20|1 . 
In ref. jHI it is proved that, in the free theory, Sj, Sj and Tj all have ranges which 
are 0(1) in units of the coarse-lattice spacing. We expect the same to hold in the in- 
teracting case. Because the original kernel Dq^{x, y) is long-ranged in physical units, 
the only way for Eq. to hold is if D~^{x,y) is long-ranged too. 

The physical interpretation of Eq. ()20|) is not completely straightforward because 
Do ^ and D-^ are not gauge invariant all by themselves. Let us examine a physi- 
cal observable. We will consider the gauge-invariant two-point function of the true 
Goldstone-boson field 7r{x). Its existence is implied by the symmetries of the one- 
component formalism ^3], and so (vr(x) vr(|/)) decays like a power of \x — y\ in the 
massless limit. We may construct this two-point function either from Dq^ or, us- 
ing Eq. (PUI) . from D~^. Had any cutoff-scale masses been induced as in the case of 
ref. jT7], would decay exponentially with a cutoff-scale rate, making it impossible 
to reproduce the correct long-distance behavior. We conclude that, whether or not 
the RG-blocked theory enjoys manifest shift and (full) hypercubic symmetries, the 
physical consequences of all the original symmetries remain intact! 

Suppose that, instead, one had started with an interacting theory obtained by 
simply gauging the flavor representation using link variables that reside on the coarse 
lattice, as in ref. J7j. Assume also that the bare quark mass is set to zero. What 
would go wrong? Denoting the new interacting Dirac operator by T>o and the lattice 
spacing by qq, The two-point function of the would-be Goldstone boson is given 
exphcitly by 

{n{x)n{y)) = {G{x,y)) , (23a) 
g{x,y) = tT{{^,^n)V,\x,y){-f,®n)V,\y,x)) . (23b) 

Since (75 ® r5)Do(75 ® ^5) = 1^1 one has Q{x,y) > 0, i.e. this correlator is strictly 
positive. The decay rate of {tt{x) vr(?/)) is now O(l/ao), because the explicit one-loop 
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calculation of ref. ^Tj shows that the quarks acquire (Lorentz-breaking) cutoff masses 
in this theory. Since Q{x, y) > 0, it is ruled out that the short range could be generated 
by destructive interference between different gauge-field configurations; rather, given 
an ensemble of configurations, we must have || T)q^(x, P ~ exp(|x — yl/ao) on each 
and every one of them. 

One could have applied RG transformations, and Eq. ()20|1 would again say that 
any long-distance physics of Vq must be reproduced by the resulting P„. However, 
we have just seen that Vq has no long-distance physics whatsoever. By Eq. (j^ . any 
number of RG transformations would not have fixed it. 

Coming back to the interacting one-component staggered theory, the algebraic 
structure of Q{x,y) is similar to Eq. ()23b|) . and again Q{x,y) > 0. In order to 
obtain the correct power-law behavior of the Goldstone-boson correlator (vr(x) vr(?/)), 
there must exist a finite probability to encounter configurations where G{x, y) has the 
same long-distance behavior. Finally, by Eq. ()20j). the long-distance behavior must 
be sustained if we use instead of D^^ to construct Q{x,y). The long-distance 
physics contained in Dq^ is faithfully reproduced by the RG-blocked D~^. 

It will be interesting to test the resulting operators numerically, and see whether 
the properties established in the free theory persist. For the RG-blocked operator Dn, 
the crucial properties are locality, suppression of the fiavor-mixing (or taste-mixing) 
part with n, and convergence to a GW operator in the massless case. As for G„ 
(Eq. (|7aj)), the question is whether this operator indeed describes only excitations 
with masses of the order of the cutoff, so that the effective action S'"g = log(det(G„)) 
is local. Last, it should be verified that the mechanism that protects the (physical 
consequences of the) symmetries of the original one-component formulation is indeed 
operative. 

On the theoretical side, an interesting idea is to make use of the notion of admis- 
sibility condition. The original concept introduced in ref. 18^ asserts that the lattice 
gauge field is constrained such that, for every plaquette, Retr(l — U^y^x)) < eo, 
where U^y{x) is the ordered product of link variables around the given plaquette, 
and eo > is a fixed (small) number. It is believed that an admissibility condition 
does not change the universality class. The utility of an admissibility condition is the 
following. Given a free lattice operator whose spectrum satisfies a certain bound, one 
expects that this bound will be modified only by 0(eo) if we promote the operator 
to a covariant one, while allowing only for gauge fields that satisfy the admissibility 
condition. 

In the case at hand, we may envisage imposing an admissibility condition on the 
link variables of the original fine lattice. This should imply that, apart from 0(eo) 
modifications, the operator bounds established in ref. will continue to hold during 
the first few blocking steps, and the same should follow for the locality properties. 

Unlike in ref. however, we now face the following problem. When the num- 
ber of blocking transformations becomes of the order of n ~ l/eo? we no longer 
have any useful bound on the gauge field. We will thus propose a stronger no- 
tion of admissibility. Considering the product of (covariant) blocking kernels Qn = 
. . . qWq(^\ we will constrain the gauge field on the original fine lattice 
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by demanding that every Wilson loop W occurring in the product QnQn satisfy 
the constraint |VV| < eo- The reasoning behind this new admissibility condition is 
the following. After n blocking steps, the length of the loops contained in QnQn will 
be 0(2'*) in units of the original lattice spacing oq. However, as explained earlier, we 
are really interested in the "two cutoff" situation where Oq = 2~"a. When measured 
in units of the coarse-lattice spacing a, all these Wilson loops have length smaller 
than some Iq = 0(1). Our new admissibility condition is therefore a natural general- 
ization of the same concept to a "two cutoff" situation. With this new admissibihty 
condition, it is plausible that the free-theory bounds will continue to hold for arbi- 
trarily large n, up to 0(eo) modifications. In other words, the deviations from the 
free-theory bounds will be 0(eo) independently of n. In addition, for very large n, the 
flavor-breaking part of the RG-blocked operator will become very small. It should be 
possible to find large enough, but finite, n, such that full fiavor ( "taste" ) symmetry 
is recovered to any desired accuracy.^ 

In conclusion, using the machinery of RG block transformations I have shown that 
the fourth-root trick is consistent for free staggered fermions. In the limit of infinitely 
many RG steps, the free staggered-fermion determinant is equal to the fourth power 
of the determinant of a one-fiavor local operator which, in the massless case, satisfies 
the GW relation, times the determinant of a local operator whose excitations have 
cutoff masses. The fourth root of the latter operator is local too. While a similar 
result for interacting staggered fermions is unlikely to be established anytime soon, 
I have discussed how to construct the fiavor representation in the interacting theory, 
while in effect maintaining all the symmetries of the one-component formalism. I have 
also suggested avenues for numerical tests, as well as a theoretical framework which 
appears to be best suited for generalizing some of the rigorous free-theory results to 
the interacting case. 
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* Since in this framework the independent hnk variables always reside on the original fine lattice, 
the limit n — > oo cannot be taken independently of the continuum limit. 
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